#include "cuda_runtime.h"
#include "device_launch_parameters.h"
#include "math.cuh"

__global__ void linearSolve(int b, float *x, float *x0, unsigned int size, float a, float c) {
    if (a == 0 && c == 1) {
		int i = (blockIdx.x * blockDim.x) + threadIdx.x;
		int j = (blockIdx.y * blockDim.y) + threadIdx.y;

		x[i*size + j] = x0[i*size + j];
		
    }
    //set_bnd(b, x);
}